Spontaneous Symmetry Breaking and Decoherence in Superconductors 
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We show that superconductors have a thin spectrum associated with spontaneous symmetry 
breaking similar to that of antiferromagnets, while still being in full agreement with Elitzur's the- 
orem, which forbids the spontaneous breaking of local (gauge) symmetries. This thin spectrum in 
the superconductors consists of in-gap states that are associated with the spontaneous breaking of a 
global phase symmetry. In qubits based on mesoscopic superconducting devices, the presence of the 
thin spectrum implies a maximum coherence time which is proportional to the number of Cooper 
pairs in the device. Here we present the detailed calculations leading up to these results and discuss 
the relation between spontaneous symmetry breaking in superconductors and the Meissner effect, 
the Anderson-Higgs mechanism and the Josephson effect. Whereas for the Meissner effect a sym- 
metry breaking of the phase of the superconductor is not required, it is essential for the Josephson 
effect. 
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I. INTRODUCTION 

Recently we have shown that spontaneous symmetry 
breaking imposes a fundamental limit to the time that a 
large spin system can stay quantum coherent. This co- 
herence timescale is t spon ~ 2nNh/ (fc^T), given in terms 
of the number of microscopic degrees of freedom N, tem- 
perature T and the constants of Planck (ft) and Boltz- 
mann (/cb).— i2j2j£j5i£j! The timescale t spon is expected to 
be a universal timescale associated with all forms of spon- 
taneous symmetry breaking, since it does not depend on 
any of the model parameters that were needed to derive 
it. In this paper we will show in detail how superconduct- 
ing systems spontaneously break their phase symmetry 
and that they have a thin spectrum associated with this^ 
The thin spectrum subtly influences the dynamics of the 
superconductor when it is used as a qubit. The resulting 
maximum coherence time is again given by the universal 
expression t spon . 

It has been proven already three decades ago by Elitzur 
that local (gauge) symmetries cannot be broken sponta- 
neously without invoking an explicitly asymmetric gauge 
fix4 Also, it has been argued recently that because the lo- 
cal gauge symmetry in superconductors cannot be broken 
spontaneously, the order should be of a purely topological 
nature, and that the low energy properties of the super- 
conducting state are determined solely by its topological 
structured At first sight then, the claim that a supercon- 
ductor possesses states related to spontaneous symmetry 
breaking that are at very low energy and within the su- 
perconducting energy gap, might come as a surprise. 

However, it is well known that the superconducting 
ground state is characterized by a definite phase and a 
corresponding uncertainty in the number of Cooper pairs. 
For a piece of superconducting material the realization of 
one specific phase is but one choice out of a manifold of 
equivalent possibilities. Its phase having just one spe- 
cific value therefore has to come about by spontaneous 
symmetry breaking. This symmetry breaking in the ther- 
modynamic limit requires the existence of a so called thin 



spectrum of total phase states whose energies all collapse 
onto the groundstate energy in the thermodynamic limit. 
Such a spontaneous symmetry breaking is not at variance 
with Elitzur's theorem because the symmetry that is bro- 
ken in a superconductor is a global U(l) phase symmetry. 
The resulting superconducting state is still manifestly in- 
variant under local gauge transformations. 

To clearly illustrate these points we will first discuss 
the superconducting state of an array of Josephson junc- 
tions. In this array the non-commutativity of num- 
ber and phase variables straightforwardly gives rise to 
spontaneous symmetry breaking and to a thin spectrum. 
Symmetry breaking in this superconductor turns out to 
be exactly analogous to the case of quantum crystals or 
magnetic systems^ After that we will switch to a mi- 
croscopic strong coupling model of superconductivity in 
which the role of gauge symmetry can be more clearly 
discussed. We will then use this model to describe a 
Cooper-pair box qubit and show that the presence of the 
thin spectrum leads to a maximum coherence time t spon 
of the qubit, which is of the order of milliseconds. Finally 
we will show how the description of the thin spectrum can 
also be incorporated into the familiar BCS description of 
superconductivity and comment on the application to 
different types of qubits. 



II. JOSEPHSON JUNCTION ARRAY 



It is well known that an array of superconducting is- 
lands coupled together by Josephson junctions can un- 
dergo a (quantum) phase transition from an insulating 
state to a superconducting statei 10 i 1:L ' 12 i 13 ' 14 i 15 The de- 
scription of a superconductor as an array of Josephson 
junctions is particularly useful to us here because it nat- 
urally focuses the attention on the the conjugate variables 
number and phase. The Hamiltonian for the Josephson 
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junction array is given by 
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Here 6j represents the phase of the superconducting or- 
derparameter of superconducting island j, while rij gives 
the number of Cooper pairs above average, and 5 con- 
nects neighboring sites. The charge or number operator 
rij = —i(d/dOj) is the variable conjugate to the phase, 
and can be written in terms of the voltage V and the ca- 
pacitance C of the Josephson junctions as rij — (C/2e)Vj. 
The coupling constants are the charging energy Ec and 
the Josephson coupling energy Ej, which are chosen to 
lie well within the superconducting regime. 

The phase 9 in this description can be thought 
of as the phase of the Ginzburg-Landau wavefunc- 
tion for the superconducting island, or equivalently as 
the phase describing the perfectly ordered BCS state 

|0> = life (kl + k|e ie 4cl fe ) |vac). 16 This phase is 
not measurable as such, but a difference in phase 
across a Josephson junction causes a supercurrent J = 
Jc sin (9 3 ■— Oj+g), and therefore phase differences are 
measurable. The condition of measurability implies the 
gauge independence of these quantities, because a gauge 
transformation by definition cannot alter the outcome of 
any experiment. The total phase is both unmeasurable 
and a gauge dependent quantity. 

The thin spectrum of the Josephson junction array 
consists of the infinite wavelength part of the Hamilto- 
nian Eq.{T]), because exactly at k = the Bogoliubov 
transformation that would diagonalize the Hamiltonian 
turns out to be singular^ This zero wavenumber part of 
H JJ which describes the collective behavior of the sys- 
tem as a whole, is given by 
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where N is the total number of superconducting islands, 
and n tot = rij is the charge of the total network of 
Josephson junctions. To see how the array can sponta- 
neously break its total phase symmetry we should add 
a symmetry breaking field to the collective Hamiltonian. 
We cannot simply add a term which involves the bare 
total phase 9 to t, because that total phase is not a gauge 
independent, measurable quantity. Instead we can look 
at the difference of phase between the Josephson junc- 
tion array and some given reference superconductor. In 
the end we will let the strength of the symmetry breaking 
field go to zero, or equivalently move the reference super- 
conductor away to infinity. The Hamiltonian including 
the symmetry breaking field thus becomes 
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For small values of A9 fot = 9 tot — 9 re f we can expand 
the cosine to quadratic order and then the Hamiltonian 



reduces to a harmonic oscillator with well known solu- 
tions in terms of Hermite polynomials, in exact analogy 
to the case of spontaneous symmetry breaking in quan- 
tum crystals and antiferromagnets. Using these Hermite 
polynomials, it is easy to show that indeed the Josephson 
junction array can spontaneously break the rotational 
symmetry of its total phase by looking at its fluctuations 
in the limit of disappearing symmetry breaking field and 
infinite number of superconducting islands: 
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Clearly the fluctuations in the total phase disappear in 
the thermodynamic limit even if only an infinitesimal 
symmetry breaking field is present. 

The symmetry broken state that is formed in that limit 
has a well defined total phase, and must thus be in a su- 
perposition of many different total number states. These 
total number states were precisely the eigenstates of the 
collective Hamiltonian Eq.f2]), which we identified as be- 
ing the thin spectrum of the Josephson junction array. 
The symmetry broken Hamiltonian also has a tower of 
low lying states that form a sort of dual thin spectrum 
which consists of all the total phase states necessary to 
build a state with a fixed total number of Cooper pairs. 
Notice that the thin spectrum states must be observable 
states, because the description of the collective dynamics 
in Hamiltonian Eq.Q is still manifestly gauge invariant. 
This also implies that the symmetry breaking which we 
have just described is not the breaking of a local gauge 
symmetry. Only the U(l) symmetry of the global to- 
tal phase is spontaneously broken, and even then only in 
the sense that its fluctuations disappear in the thermo- 
dynamic limit, so that its value relative to that of some 
other, external superconductor will be fixed. 

The fact that we needed to introduce an external su- 
perconductor as a deus ex machina to fix the phase of 
our Josephson junction array should come as no surprise. 
The situation is in fact precisely analogous to that of 
breaking the translational symmetry of a crystal. In that 
case one can only assign a definite value to the position of 
the symmetry-broken crystal by measuring the distance 
of its center of mass to some external reference point 
(which in that case can the observer himself). The posi- 
tion of the entire system of crystal and observer together 
is still completely arbitrary (or at least unmeasurable for 
the observer), even in the symmetry-broken state. 



III. LOCAL PAIRING SUPERCONDUCTOR 

From the previous section it is clear that that the non- 
commutativity of number and phase naturally gives rise 
to the presence of a thin spectrum in a superconducting 
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system. Now we will examine how the superconducting 
orderparameter comes about by spontaneous symmetry 
breaking in the first place and how this is related to gauge 
symmetry. This relation was not visible in the context 
of a Josephson junction array, because there we started 
out with islands that were postulated to be in a super- 
conducting state. That way we could describe the whole 
system with an effective Hamiltonian that only consisted 
of observables related to the macroscopic properties of 
the superconducting state. 

For a more general description of superconductivity 
we start out with a microscopic Hamiltonian for a single 
superconductor that incorporates the effects of the gauge 
field. The simplest such model is the extensively studied 
local pairing, negative-!/ Hubbard mode l 17 ! 18 
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The z projection of the pseudospin measures the local 
electron density, while the xy components provide the dy- 
namics of the Cooper pairs. Adding a chemical potential 
fi that determines the overall electron density and thus 
explicitly breaks the electron-hole symmetry, we find the 
effective Hamiltonian 
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Here J is defined to be 2t 2 /\U\, and h = \U\ — 2fx deter- 
mines the overall electron density. Away from half filling 
the global SU(2) symmetry of the Hamiltonian is man- 
ifestly broken, and what remains is the U(l) symmetry 
that describes rotations around the z-axis. It is the spon- 
taneous breaking of this global U(l) symmetry that will 
yield the superconducting state. 



Here c] creates an electron on site j, S connects neigh- 
boring sites and rij counts the number of electrons. The 
reason to consider this local pairing model rather than for 
example the BCS model for superconductivity is the fact 
that this model is explicitly gauge invariant, while the 
BCS model is not. From the symmetry point of view, 
the models are the same: there is no phase transition 
in going from weak to strong coupling superconductivity, 
only a cross-over. If we parametrize the hopping in terms 
of a uniform amplitude and a bond dependent phase as 
tj = te^i , then minimal coupling allows us to identify 
the phase of the hopping parameter with the electromag- 
netic vector potential integrated along the bond under 
consideration, so that ^ = ^ Jj +S A s \t)dt. Thus the 
Hamiltonian is invariant under the gauge transformation 



i-tfU)J 
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which immediately implies 
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We focus on the strong coupling limit where U t, 
so that we only need to consider the physics of the lower 
Hubbard sector. On each site there will thus be either 
a pair of electrons or no electrons at all. Single elec- 
tron excitations are only virtually allowed and give rise 
to pair-pair interactions. The effective low energy Hamil- 
tonian is given by a second order perturbation expansion 
in the hopping and can be written in terms of pseudospin 
operators that are defined by 
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IV. MEISSNER EFFECT AND 
ANDERSON-HIGGS MECHANISM 

Before we discuss the actual spontaneous symmetry 
breaking and the thin spectrum associated with it, we 
will show that already on the semi-classical level, the 
model Eq.@ can be seen to expel magnetic field lines 
from its ground state and to give propagating electro- 
magnetic modes in its bulk a finite effective mass (the 
Meissner effect and Anderson- Higgs mechanism). 

To find a semiclassical description for the groundstate 
of the S = 1/2 pseudospin Hamiltonian H e f f , we intro- 
duce generalized coherent states of the form 



\^ class) = ]J ( e ^ Sm ( ^ 
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In this expression the angles an d Oj are the Eulcr 
angles which describe the classical vectors that replace 
the quantum spins in the classical state. To find the 
semiclassical groundstate energy we need to minimize 
the expectation value of H e ff in the generalized coher- 
ent state with respect to the orientations of the classical 
spin-vectors. It is easy to check that the classical energy 
can be minimized by first fixing the azimuthal angles 9j 
uniformly throughout the system, after which the energy 
up to a constant is given by 
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where p is a constant set by the optimized value of 6. In 
this expression the global SU(2) rotational symmetry has 
already been broken explicitly by the effect of the field h 
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(which fixed the azimuthal angles) . What is left is the in- 
plane U(l) rotational symmetry of the polar angles 4>j- 
The classical state with lowest energy links the values 
of these polar angles to the bond variables ipj. These 
variables are in turn connected to the electromagnetic 
vector potential. One finds that the condition for the 
angles <j>j that minimizes the energy is 
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where Aj is the average value of the vector potential 
along the bond. At distances much larger than the lattice 
spacing a this expression becomes 



A(r) = gv0(r) 



(13) 



The classical groundstate is thus a state in which the 
electromagnetic potential is proportional to the gradient 
of the scalar field <p, which of course immediately im- 
plies that the rotation of A will vanish, and thus that 
the condensate does not allow any magnetic field to pen- 
etrate its bulk -a clear indication that the semiclassical 
groundstate is indeed a superconducting state. 

For a full description of the Meissner effect, however, 
it is not enough to show that the semiclassical ground 
state does not contain a magnetic field. One also needs 
to demonstrate that the superconducting classical state 
will also actively expel externally applied electromagnetic 
fields; i.e. that the electromagnetic excitations in the 
system are gapped and massive. To do so we add an 
external electromagnetic field to the semiclassical energy 
density: 
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This expression can be simplified by introducing a new 
vector field A' = A— V</>. Because the newly defined field 
A' is formally equivalent to a gauge-deformed version 
of the electromagnetic field A, we can be sure that the 
electromagnetic energy E 2 + B 2 looks the same in terms 
of A' as it did in terms of A. If we take the continuum 
limit and expand the cosine to second order, we thus find 



E class ~ J' (A') 2 + ~ (V x A') 2 + ~ ( A' 



(15) 



Note that both expressions for the energy given above 
are fully gauge invariant. Using the Hamilton equations 
it can immediately be checked that the latter expression 
for the classical energy yields only massive propagating 
modes in terms of the field A'. Due to this Anderson- 
Higgs mechanism, the physical excitations of the system 
-which are combined modes of the electromagnetic field 
and phase degree of freedom- have a finite energy gap. 
This prevents an external electromagnetic field to pene- 
trate the bulk of the superconductor. 



V. THIN SPECTRUM 

We have seen in the previous section that the sem- 
icalssical groundstate of the tight binding negative-C/ 
Hubbard model displays the Meissner effect and gen- 
erates a mass for the electromagnetic modes via the 
Anderson-Higgs mechanism. In fact in that semiclassical 
description the question as to whether the global U(l) 
symmetry was broken or not never arose. We thus con- 
clude that the Meissner effect and Anderson-Higgs mech- 
anism occurs regardless of whether the superconductor 
has a well-defined total phase. 20 ' 21 ! 22 The situation is 
similar to that in antiferromagnets, where long range an- 
tiferromagnetic correlations exist both in the symmetric 
singlet ground state and in the symmetry broken Neel 
state. 

To see the effects of spontaneous symmetry breaking 
in the negative-C/ Hubbard model we need to describe 
the formation of the symmetry broken ground state in 
a more analytical manner by studying the exact eigen- 
states of the collective part of the Hamiltonian without 
resorting to semi-classics, just as we did for the Josephson 
junction array. The difficulty in such a global description 
will be to correctly account for the gauge field, which can 
fluctuate locally. To circumvent this problem we intro- 
duce transformed pseudospins, analogous to what is done 
in the weak coupling theor3 r 18 i 19 
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The summation in the exponent is over some path con- 
necting position j to some origin j = 0. For simplicity 
we will assume the applied external magnetic field to be 
zero from here on. Notice that the individual transformed 
pseudospin operators of equation Eq. (|16p are not gauge 
invariant. Their purpose is to transform the local gauge 
transformations of the actual pseudospins S into a global 
transformation of the new pseudospins a: 

a + e - 2i ^^oM' + MfO'+S)-f(3'))] e 2i^J(j) s + 



(17) 



A gauge transformation therefore corresponds to a ro- 
tation along the z-axis of all pseudospins by the same 
angle. It is this global character of the gauge transfor- 
mations on the transformed pseudospins that allows us 
to switch to a description of just the collective behavior 
of the system without invoking any specific gauge choice. 
In terms of the transformed pseudospins the effective low 
energy Hamiltonian Eq.@ becomes 



H eff = J Y, a 3 
3,6 



7 3+s- h J2 a l 



(18) 



The parameter J is positive because the electron pairs re- 
pel each other. This Hamiltonian therefore describes an 
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FIG. 1: A schematic one-dimensional representation of the 
classically realized state of the tight binding superconductor. 
The arrows are a classical cartoon for the transformed pseu- 
dospins a. 



antiferromagnetic interaction between neighboring pseu- 
dospins, and an overall, uniform magnetic field. The clas- 
sical state that we expect to find in terms of the pseu- 
dospins a is therefore a canted antiferromagnet. That 
is, an antiferromagnet in which all spins are uniformly 
canted out of the z = plane, but in which the xy- 
projections still form an antiferromagnetic pattern (see 
figure [1]). As we noticed before, the canting of the spins 
which breaks the full SU(2) symmetry down to U(l) is 
done explicitly by the field h, while the breaking of the in- 
plane U(l) symmetry into an antiferromagnetic structure 
will have to be realized through spontaneous symmetry 
breaking. 

The thin spectrum of the Hamiltonian H e f f consists 
of the states necessary to construct the symmetry bro- 
ken classical state. These thin spectrum states describe 
the dynamics of the superconductor as a whole, and as 
before they can be found at the singular points of the Bo- 
goliubov transformation that diagonalizes the quadratic 
part of the Hamiltonian^ In the antiferromagnet both 
the point k = and the point k = tt are singular. The 
resulting collective part of the Hamiltonian is given by 



4J 

77 



- hai 



(19) 



where <Ja,b denotes all spins on the A, B sublattice and 
a tot is the combination of all spins on the entire lattice. 
The modes that form the thin spectrum are not affected 
by the Anderson-Higgs mechanism that we discussed be- 
fore because they are zero- wavelength, global excitations 
which only affect the system as a whole. The coupling 
of the phase degrees of freedom and the electromagnetic 
field in the Anderson-Higgs description exists only at fi- 
nite wavelength, as can be easily seen in equation Eq. (fT^|) 
because a global transformation fa — > fa + 5(f) for all i 
leaves the coupling of the phase to the electromagnetic 
field invariant. The Hamiltonian H co u therefore captures 
all the relevant collective, low energy behavior of the 
model. 

The collective Hamiltonian Eq. (fl9| is just a Lieb Mat- 
tis model in a uniform magnetic field,— i 23 i 24 i 25 and the 
eigenstates are readily identified as the states labeled by 
the quantum numbers a A ,a B , a tot and of ot . The differ- 
ence between this collective model and the one describing 
the spontaneous symmetry breaking in antiferromagnets 
is the field h, which reduces the symmetry from SU(2) 
to U(l). The ground state has maximum total spin on 
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FIG. 2: Schematic representation of the dispersion relation of 
the low-energy, low-momentum states of a finite size super- 
conductor: the states at finite k are gapped, while the thin 
spectrum states at k — lie within the gap. 



both the A and B sublattice, and has a total spin that is 
equal to its z-projection of af ot = a tot = (hN)/(4J). Ex- 
citations of the quantum numbers a A and a B are gapped 
with an energy J from the groundstate, because of the 
infinitely long range of the interactions in H co u . We will 
henceforth set these quantum numbers to their maximum 
value and only consider the low energy excitations which 
describe the behavior of the entire system as a whole. We 
can relabel the eigenstates by introducing 
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Here a is the groundstate value for the z-projection of the 
pseudospin: a — {hN)/{AJ). In terms of the quantum 
numbers n and y, the effective Schrodinger's equation 
for the excitations of the Lieb-Mattis Hamiltonian H co u 
becomes 
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From this equation it is clear that the excitations labeled 
by n will play the role of the thin spectrum for the tight 
binding superconductor. It can be easily checked that 
indeed the contribution of these states to the partition 
function vanishes in the thermodynamic limit. The exci- 
tation labeled by y on the other hand is a collective exci- 
tation that changes the z-projections of all pseudospins 
and costs an energy proportional to the chemical poten- 
tial to excite. This is, in other words, the quantum num- 
ber that determines the average total number of Cooper 
pairs in the superconductor. 



VI. BREAKING THE SYMMETRY 

To study the spontaneous symmetry breaking of H co u 
we will have to introduce a symmetry breaking field, 
which we will send to zero again at the end of the calcula- 
tion. In analogy to the symmetry breaking term that we 
used in the Josephson junction array, we will again intro- 
duce a second external superconductor which is weakly 
coupled to the first, and let the coupling tend to zero: 
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Here denotes the pseudospin operators in the exter- 

nal superconductor and H ext describes its dynamics. For 
simplicity we assume that both the Hamiltonian H co u 
and H ext are given by the collective model of equa- 
tion Eq. (fT9|) . Notice that the coupling Hamiltonian is 
still completely invariant under local gauge transforma- 
tions which act on the pseudospins ga,b as well as on 
Tia,b- The symmetry that is broken in equation Eq. (122p 
is the U(l) phase symmetry that rotates all pseudospins 
a, but keeps the pseudospins X fixed: the phase differ- 
ence between the superconductors can acquire a finite 
expectation value in the symmetry broken state, without 
violating Elitzur's theorem. 

To explicitly see the influence of the external supercon- 
ductor, we would liketo use the total phase of the system 
H ex t as a reference point for measuring the phase of H co u . 
To do so, we need to temporarily abandon the manifest 
gauge invariance of our description, and assume that the 
total-phase symmetry of H ext has already been broken. 
This allows H ext to generate a symmetry breaking field 
for the cr-pseudospins subsystem. Using for this (mean) 
field the expectation values (E A ) = — (E B ) = and 
(E A ) = (Eg) = 0, the effective Hamiltonian for the orig- 
inal superconductor reduces to 



monic oscillator equation, 
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with B = TNs 7 and the number of pseudospins in 
the external superconductor. Notice that the symmetry 
breaking term in this Hamiltonian is not gauge invari- 
ant. As a gauge transformation corresponds to a uni- 
form rotation of all spins on the entire lattice around the 
z axis, any explicit choice for the direction of B along a 
particular axis in the iry-plane is connected by a gauge 
transformation to any other directions in the plane. By 
fixing the total phase of the external superconductor to 
lie along the x-axis in equation Eq. (|23|) . we have thus 
implemented a specific gauge choice, and we will have 
to check afterward if the conclusions based on calcula- 
tions in this particular gauge fix are robust under gauge 
transformations. 

The matrix elements of the symmetry breaking field 
in the basis \n, y) can be computed by performing a sum 
over Clebsh-Gordon coefficients^ In the limit (n, y) <C 
(a,N), the Hamiltonian can be written in terms of its 
matrix elements as 
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where f(y) = (2 - g) \J (f) - a 2 . If we write 
the eigenfunctions of this equation as \x, y) = 
^ ( n ' x ) l n ' V) an d take the continuum limit, then 
Schrodinger's equation reduces to the well known har- 
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With ^ = BNf(y) aild V 
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The 



wavefunctions *S? are the eigenfunctions of the harmonic 
oscillator, which can be written explicitly in terms of Her- 
mite polynomials. The corresponding eigenvalues obey 
v = (x + l/2)w, and thus we find the energies of the sym- 
metry broken collective Hamiltonian Eq. (|2"3"|) to be given 
by 



E(x,y) = E° coll 



hy - ^BNg(y) 



+ (x + - V2JB y/tfy), (26) 



where g(y) = (l - 2f ) yjl - (jj) ■ The term cx BN 
in this expression shows that the symmetry of the sys- 
tem will be spontaneously broken: even if only an in- 
finitesimally small symmetry breaking field is present, 
the pseudospins can gain an infinite amount of energy 
in the limit of N — > oo by aligning with that field. In 
the thermodynamic limit the alignment will thus happen 
spontaneously and the resulting symmetry broken state 
is exactly the expected canted antiferromagnet. 



VII. GAUGE VOLUME 

Having found the the eigenfunctions of the collective 
symmetry broken Hamiltonian, the question arises what 
these states represent, and even if they are truly physi- 
cal states. As mentioned before, the symmetry breaking 
field in the collective Hamiltonian Eq. (|23|) acts as an im- 
plicit gauge fix. It is not a priori clear whether or not 
this (non-physical) gauge fixing introduced any extra un- 
physical states in the spectrum. If we define the gauge 
volume of a certain state to be the collection of all states 
that are connected to it by a gauge transformation, then 
making a specific gauge choice in the Hamiltonian can 
in principle lead to the erroneous identification of states 
within the same gauge volume as seperate physical states. 
The question is thus whether the excited states of 
that we found are part of its ground state gauge volume 
or not. 

The ground state of the collective Hamiltonian is an 
ordered antiferromagnet in terms of pseudospins, and we 
have seen that it corresponds to a superconducting state 
of Cooper pairs. The excitations labeled by x in the 
pseudospin picture must involve the superposition of col- 
lective excitations with wavenumbers k = and k = ir. 
However as mentioned before, the gauge volume of this 
system is made up of global uniform rotations of the en- 
tire pseudospin lattice around the z-axis. Proving that 
these excitations are not within the gauge volume of the 
groundstate wavefunction, therefore amounts to showing 
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X=10 




FIG. 3: The overlap between the thin spectrum state \x) and 
the rotated groundstate R (9) |0), as a function of the angle 
of rotation 8, for different values of x. To make this graph we 
used the values J = 10, B = h = 1 and N = 100. 



that the excited states cannot be written as only a global 
rotation of the groundstate. Using the explicit formulas 
for the eigcnfunctions of it is easy to check that in- 
deed the overlap between the state \x = X) and the state 
| a: = 0), rotated over an angle 9, is one if and only if both 
X and 9 are zero (see figure \S§ . This proves that indeed 
the excited state cannot be written only as a global rota- 
tion of the groundstate, and thus that the excited state 
is not within the ground state's gauge volume.— 



VIII. JOSEPHSON EFFECT 

We have seen that an infinitesimally weak coupling be- 
tween the local pairing superconductor and a second ex- 
ternal superconductor gives rise to spontaneous symme- 
try breaking of the total phase difference between the two 
superconductors. In this section we show that this sym- 
metry breaking is an essential prerequisite for the obser- 
vation of the Josephson effect. This effect is well-known 
to rely on the existence of a phase difference between 
two superconductors, in which case a finite (but weak) 
coupling between the two systems causes a supercurrent 
between themi 26 i 27 i 28 

To study the Josephson effect in the local pairing 
model for superconductivity, we will follow the approach 
of Anderson by treating the hopping term in equa- 
tion Eq. (|22[) as a small perturbation and computing the 
change in energy due to iti2S By taking the derivative of 
the the first order correction to the energy with respect 
to the vector potential we can then directly obtain the 
supercurrent. 

First we construct the state in which there are no ex- 



citations y in either superconductor, and in which the 
thin spectrum state of both superconductors is \x = 0), 
which we denote as |0) = |0, 0) s |0,0) ff . Now we rotate 
the phase of the S-superconductor by an angle 9 and 
get the state \9) = R(6) |0), where R rotates the state 
|0,0) s over an angle 9 around the z-axis. The fist order 
correction to the energy is 



= T 



T(0|(S+^+E+^) 
N a Nx 



H.c. 



■ cos (9) 



* 0' + 1, 0) * (n - 1, 0) * (n', 0) * (n, 0) 

n,n' 

= AEcos(9), (27) 

The summation over the Hermite polynomial wavefunc- 
tions \& (n, x) can easily be evaluated numerically and is 
of order unity. In the last line the energy change AE 
is implicitly defined. The expression above for the en- 
ergy difference in terms of the phases of the transformed 
pseudospins a and S is gauge invariant, because a gauge 
transformation rotates all pseudospins simultaneously, 
and leaves the phase difference 9 invariant. As we wish 
to determine the derivative of AEi with respect to the 
vector potential, we need to go back to the formulation 
of the problem in terms of the original S pseudospins. By 
inspection of Eq. ([27]) it is clear that the gauge-dependent 
expression for the energy change will be 



AE 1 = AEcos[9- — 

he 



A ■ dl 



(28) 



where 9 describes the phase difference in terms of the 
original pseudospins, and the integration runs over a line 
connecting the two superconductors. We can now di- 
rectly find the supercurrent by taking the derivative of 
the total energy with respect to the vector potential. Do- 
ing so and then transforming back to the pseudospins a 
and E yields 



5(H) 
' SA 



2eAE 



sin (0) , 



(29) 



which is precisely the expected Josephson current. 28 No- 
tice that the occurrence of the Josephson effect is a direct 
consequence of the phase symmetry breaking. Starting 
from a symmetric, total-number ground state for the su- 
perconductors, the first order correction to the energy 
induced by the hopping term vanishes, so that there is 
no Josephson current in that case. 

This observation raises the question what happens in 
practice when we try to bring two isolated pieces of su- 
perconductor closer together and allow them to weakly 
couple. Assuming that the superconductors were so well 
isolated that their initial states were total number states, 
the first influence of the coupling will be to break the 
phase symmetry and cause the superconductors to ac- 
quire a well-defined phase difference. If two supercon- 
ductors approach each other from infinity the actual 
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value of the phase difference 9 will be zero -due to the 
energy gain Eq.(|2"T|)- and no Josephson current of the 
type Eq. ([29|) is present. However, a finite phase differ- 
ence between the two pieces of superconductor can be 
induced by first applying a voltage difference between 
the two systems and then letting them approach. In the 
Hamiltonian the voltage bias can be introduced into the 
hopping term by the Peierls construction: 



jo 



A w A 



^B a B 



H.c. 



(30) 



where 8 depends on the applied voltage. This form of the 
coupling term causes the symmetry to be broken such the 
phase difference between the superconductors is 6. If the 
bias is switch off again, a Josephson current as in Eq. J29J) 
is induced. 



IX. DECOHERENCE 

We would now like to apply the results of the previ- 
ous sections to the description of quantum coherence. In 
analogy to the result for antiferromagnets, 1 we expect 
the existence of the unobservable thin spectrum to give 
rise a maximum coherence time t spon oc Nh/ksT. 

Let us define a qubit made of the eigenstates of the 
collective part of the local pairing superconductor. If 
temperature is sufficiently low (i.e. ksT <C J, h) then 
we can use the states y = and y = 1 as the compu- 
tational states of such a qubit. These states correspond 
to states with a different number of Cooper pairs, and 
qubits of this type have been made experimentally in the 
form of Cooper-pair boxes. 29 ' 30 ! 31 In these Cooper-pair 
boxes a superconducting island can be brought into a su- 
perposition of having N and N + I Cooper-pairs present. 
Superpositions of this type can reach coherence times of 
up to 500 ns-2 2 .^ 

In our local pairing description of the qubit, the initial 
state of the system must be a thermal mixture of thin 
spectrum states. After all, controlling these states ex- 
perimentally is practically impossible. The initial state 
should then be brought into some superposition of the 
computational states y = and y = 1, so that it can be 
used in a quantum computation. Because we know all 
eigenstates and eigenvalues of the Hamiltonian exactly, 
we can then explicitly follow the time evolution of the 
superposition^ The complete process is thus described 
by the time dependent density matrix 

Pt<o = ^e-W)M><*,0| 

p t =0 = ^^e-^ o )[|x,0) + | 2 :,l)][(a ; ,0| + (x,I|] 

X 

Pt>o = ^^2e^ E ^\\x,0)(x,Q\ + \x,l)(x,l\ 

X 

+e -|(B(x,0)-E(x,l))t ^ J] + h ( , 1 _ (31) 



where Z is the partition function at t < 0. The thin 
spectrum states labeled with x cannot be observed or 
controlled experimentally, and they should therefore be 
traced out of the final density matrix. The remaining 
reduced density matrix then shows the coherence of only 
the superposition of y states. The disappearance of the 
off-diagonal matrix element of the reduced density matrix 
serves as a measure of the resulting coherence time, and 
it can easily be checked that this coherence time is given 
by 



t. 



2ttH a 



(32) 



Here a signifies, as before, the average number of Cooper 
pairs on the superconducting island in the groundstate. 
This coherence time is the maximum coherence time of a 
superconducting island, which is limited by the existence 
of a thin spectrum in the superconductor. Just as in the 
cases of crystals and antiferromagnets, the details of the 
model (e.g. J or h) do not enter into the expression for 
the maximum coherence time, which thus looks like a 
universal timescaleji^ 

Filling in the values for the constants h and ks and 
taking a ~ 10 6 and T ~ 40 mKp- we find a maximum 
coherence time for the experimentally realized Cooper 
pair boxes of ~ 0.5 ms. Clearly this timescale set by the 
presence of the thin spectrum states which are associated 
with the spontaneous symmetry breaking, is much larger 
than the timescale that is the current limit to coherence 
of the Cooper-pair boxes due to other environmental fac- 
tors. However, it is well possible that the limit set by the 
thin states will come within the experimental reach in 
the near future, either because the isolation from exter- 
nal sources of decoherence will be developed further, or 
because the size of the Cooper-pair box itself is reduced 
even more. 



X. BCS SUPERCONDUCTOR 

In the previous sections we have shown that the super- 
conductive groundstate is a state with a spontaneously 
broken U(I) symmetry. As a consequence the supercon- 
ductor must have a thin spectrum of states that describe 
the collective excitations on top of the ground state. In 
the case of a local pairing model for superconductivity 
we have found an explicit expression for these thin states 
and we have shown how they can cause decoherence if we 
try to use a superconductive island as a qubit. 

It could be argued that the local pairing model is some- 
what pathological, and not really representative for real- 
life superconductors, even though from the point of view 
of symmetry the model is equivalent to a weak coupling 
model (because there is no phase transition which sep- 
arates the two). We will therefore also work out the 
symmetry breaking and decoherence in a BCS descrip- 
tion, and show that although the details of the picture 
change, the underlying physics is exactly equivalent, and 



in fact gives rise to the exact same conclusions regarding 
the thin spectrum and the timescale on which decoher- 
ence will set in. The draw-back of doing the calculation 
in a mean field BCS description is that it cannot be done 
in a manifestly gauge invariant way, so that the role of 
the vector potential is obscured. 

After creating Cooper pairs, we arrive in the standard 
BCS theory at the effective Hamiltonian 



H 



BCS 



£ e * (< 



t. 



■ J2 u r 'l f;t // : 

k=£k' 



(33) 



Here we have adopted the convention to write (fc, f) as 
k and (— fc, |) as —k. The dispersion of the bare Fermi- 
sea is characterized by while U is the effective pairing 
interaction due to phonon exchange. U is non-zero and 
attractive only in a shell around the Fermi energy with a 
width of about the Debije energy. It is easy to see that 
extensivity of the model in fact requires U to be inversely 
proportional to the total number of electrons in the sys- 
tem. We will therefore redefine the pairing potential as 
U = V/N, where N denotes the total number of electrons 
in the k-space shell in which U is non-zero. 

By writing down the Hamiltonian Eq. (|33[) we have as- 
sumed that there is no external magnetic field and we 
have fixed the gauge to ensure that the electromagnetic 
vector potential vanishes everywhere. Anderson showed 
that the BCS Hamiltonian in this form can be rewritten 
as a spin problem by introducing the pseudospins 19 



si 



1 c k c k c -k°-k 



(34) 



In the subspace without any quasiparticles (i.e. = 
ri—k Vfc), the Hamiltonian up to an overall constant be- 
comes 

Hbcs = "2 " £ {S%St, + SlSl ) . (35) 



k^k' 



Interpreted at face value, this Hamiltonian describes 
pscudo spins on a lattice which has position-label k. On 
this lattice, three different and independent regions can 
be identified. In the region k < kp — kp, (where /cf is the 
Fermi wavenumber and kp, the Debije wavenumber) we 
know that the pairing potential vanishes and is neg- 
ative, so that all pseudospins in that region will point 
down, which corresponds to completely filled electronic 
states. In the region k > hp + kp, the pairing potential 
is zero as well, but here tk will be positive, causing all 
spins to point up, and all electronic states to be empty. 
In the shell of width kp, around kp a more interesting sit- 
uation occurs. There V is nonzero (and approximately 
constant), while switches sign right at kp. The pseudo- 
spin structure that one would classically expect in that 




FIG. 4: A schematic representation of the region of width ko 
around kp. The arrows represent the pseudospins S. Sponta- 
neous symmetry breaking causes the projections of the pseu- 
dospins in the a;y-plane to align. 



region is that of a magnetic domain wall: the pseudospins 
point up at one end of the region, then continuously fall 
over until they reach the xy plane exactly at kp, and then 
they continue on until they point down at the other end 
(see figure |3]). Electronically that structure corresponds 

to the BCS wavefunction \\ k (u^ + v^c^c^^j |vac). 

The Hamiltonian Hbcs however is invariant under ro- 
tations around the z-axis, and the exact groundstate will 
also obey this symmetry and have a completely delocal- 
ized projection of the pseudospins on the xy plane. To 
form a true domain wall, and thus the classical super- 
conducting state, this U(l) symmetry will have to be 
spontaneously broken. 

Because the symmetry breaking will only have an ef- 
fect in the region around kp and because this region is 
fully decoupled from the other two regions of fc-space, we 
will focus solely on that shell from now on, and define 
all sums over k to run from kp — kp, to kp + kp,- The 
collective dynamics of the system will again be described 
by the singular points of the Bogoliubov transformation 
which diagonalizes the Hamiltonian. Because of the fer- 
romagnetic sign, the collective model in this case consists 
of only the k = part of equation Eq. ([55]) : 



H, 



coll 



N 
V 
N 



qz 
Hot°tot 



V(N-l) 
N~ 2 



qx qx _i_ qy qy 
J tot J tot ' a tot°t, 



[St 



QZ QZ ~| 

D tot D tot\ 



(36) 



where Stot = X)fc Sk and where in the last line we have 
neglected terms of order 1/N 2 and set t t ot — 0. The lat- 
ter can be thought of as a strong coupling approximation, 
in the sense that the Hamiltonian Eq. (f36|) will certainly 
be relevant in the region where e tot <C V. We will discuss 
different approximations for e to t at the end of this section. 
The eigenstates of the collective Hamiltonian are trivially 
found to be labeled by the total spin quantum number 
S and its z-projection M, while the corresponding ener- 
gies are given by E co u(S, M) = -V/N (S {S + 1) - M 2 ) . 
The thin spectrum in this case is labeled by M, and de- 
scribes states with different total electron densities. The 
total spin excitations labeled by S on the other hand, are 
gapped with an energy ~ V. To break the xy-symmetry 
of H co u we can add a symmetry breaking field —BSf ot 
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along for example the x-axis. After evaluating its matrix 
elements^ and taking the continuum limit, Schrodinger's 
equation can once again be written as a harmonic oscil- 
lator equation 

- ( M > x ) + \u 2 MH (M, x) - i/* (M, x) ,(37) 

With w 2 = ^_ andi; = 1+ E(S,*)-jWS,0) _ sym _ 

metry broken wavefunctions \S,x) = J2 M &(M, x) \S, M) 
thus have energies 

E(S,x) = ~S{S + 1)-BS 

+ H)W¥ < 38 » 

In the ground state S will be maximal (i.e. N/2), and 
then the term oc NB in the energy signals spontaneous 
symmetry breaking: in the thermodynamic limit the sys- 
tem can gain an infinite amount of energy by aligning 
with an infinitesimally small symmetry breaking field. 
The collective excitations that make up the (dual) thin 
spectrum on top of the symmetry broken ground state are 
labeled by x. Their energies are slightly influenced by the 
remaining collective quantum number S. If we make a 
superposition of total spin states and trace away the un- 
observable thin spectrum, then this small shift in the thin 
spectrum's energy levels will cause the decoherence of the 
visible reduced density matrix, in a manner completely 
analogous to the one described in equation Eq. (j3Tj) . The 
resulting maximum coherence time is given by 

2itH 

tspon = ~, Tfi 1 (39) 

k B T 

where N counts the number of states in the k space vol- 
ume of kp> around hp, which is proportional to the num- 
ber of Cooper pairs in the superconducting condensate. 
So we find the same universal form for the expression of 
the coherence time set by spontaneous symmetry break- 
ing as in the case of the local pairing model for super- 
conductivity. 

As mentioned before, the collective Hamilto- 
nian Eq. (|36| can be seen as a strong coupling limit, 
because we require e tot to be much smaller than V. We 
can drive the system to a somewhat weaker coupling 
regime by reincluding an approximate form of J2k e kS^ 
into H co u . One possible choice for such a term would be 
t (iSf . — ), which acts as a boundary condition, 
pulling the pseudospins down at the low k boundary 
and up at the other end. A second choice could be the 
inclusion of the term t (S A — S B ) where Sa consists of 
all spins with k < kp and Sb denotes spins above the 
Fermi surface. In the latter case we should take care 
that t cannot be too great, for if it would dominate 
over V everywhere, then it would transform the domain 
wall structure of the superconducting state into a trivial 
Fermi-sphere structure again. It turns out that after 



some elaborate algebra both of the above cases give 
the exact same form for the thin spectrum and the 
maximum coherence time as the "bare" model H co u did. 



XI. CONCLUSIONS 

In this paper we have shown that the non-commuting 
observables of number and phase in a superconductor 
give rise to spontaneous symmetry breaking and an asso- 
ciated thin spectrum. We have given explicit expressions 
for these thin spectrum states in an array of Joseph- 
son junctions, in a tight-binding, negative-// Hubbard 
model, and in the BCS model for superconductivity. Us- 
ing the negative-C/ Hubbard model we have commented 
on the relation between spontaneous symmetry breaking 
and its associated thin spectrum, and the Meissner ef- 
fect, the Anderson-Higgs mechanism and the occurrence 
of Josephson currents. For the occurrence of the Meissner 
effect the phase symmetry actually need not be broken, 
but for the Josephson effect it does. We have also given 
a description of a gedanken experiment in which the su- 
perconductor is to be used as a qubit, and we have shown 
that the presence of the thin spectrum states associated 
with spontaneous symmetry breaking will lead to deco- 
herence of the qubit within the time t spon — 2TthN/kBT, 
where N counts the number of Cooper pairs involved. 
This result was obtained in the negative-?/ Hubbard 
model as well as the BCS model. The timescale t spon is 
universal in the sense that it does not depend on the un- 
derlying model parameters. Its form coincides precisely 
with that of the decoherence time induced by thin spec- 
trum dynamics in antiferromagnets and quantum crys- 
tals. 

The maximum coherence time that we found here for 
superconducting devices should apply directly to experi- 
mental realizations of the so called Cooper pair box, and 
thus give a maximum coherence time of the order of mil- 
liseconds. The decoherence caused by the thin spectrum 
at the moment is much weaker than that caused by other 
sources, but it may well come within experimental reach 
in the near future. To apply the results of this paper to 
other types of superconducting qubits, such as for exam- 
ple superconducting flux qubits, one should adjust the 
models used here in order to also accommodate for the 
presence of an external magnetic flux and an associated 
supercurrent in the groundstate. Because of its universal 
nature it is expected that the decoherence time set by 
the thin spectrum in these cases also will be given by the 
timescale t spon — 2ixtiN jksT . 
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